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We study coherent backscattering of a monochromatic laser by a dilute gas of cold two-level 
atoms in the weakly nonlinear regime. The nonlinear response of the atoms results in a modifica- 
tion of both the average field propagation (nonlinear refractive index) and the scattering events. 
Using a perturbative approach, the nonlinear effects arise from inelastic two-photon scattering pro- 
cesses. We present a detailed diagrammatic derivation of the elastic and inelastic components of the 
backscattering signal both for scalar and vectorial photons. Especially, we show that the coherent 
backscattering phenomenon originates in some cases from the interference between three different 
scattering amplitudes. This is in marked contrast with the linear regime where it is due to the 
interference between two different scattering amplitudes. In particular we show that, if elastically 
scattered photons are filtered out from the photo-detection signal, the nonlinear backscattering en- 
hancement factor exceeds the linear barrier two, consistently with a three-amplitude interference 
effect. 

PACS numbers: 42.25.Dd, 32.80-t, 42.65-k 



I. INTRODUCTION 

Propagation of light waves in disordered media is an ac- 
tive research area since hundred years ago now. The orig- 
inal scientific motivation came from astrophysical ques- 
tions about properties of light radiated by interstellar 
atmospheres [J, |2| • Then, within the first decades of the 
twentieth century, the foundations of light transport in 
this regime were laid, leading to the radiative transfer 
equations 

II II 13 The basic physical ingredient of 
these equations is a detailed analysis of energy trans- 
fers (scattering, absorption, sources, etc). Sufficiently 
far from any boundaries, the long-time and large spatial 
scale limits of these equations give rise, in the simplest 
cases, to a physically appealing diffusion equation. 

One important feature of this theory is to consider that 
any possible interference effects are washed out under dis- 
order average. This is a random-phase assumption. For 
a long time, it was believed that this was still the case 
on average for monochromatic light elastically scattered 
off an optically thick sample even if, for a given disorder 
realization, one observes a speckle pattern Q indicating 
that phase coherence is preserved by the scattering pro- 
cess. Theoretical and experimental works in electronic 
transport [H, 0, 0] made soon clear that this random- 
phase assumption was wrong in the elastic regime. De- 
pending on the disorder strength, partial (weak local- 
ization regime) or complete (strong localization regime) 
suppression of diffusive behavior has been predicted, pro- 
vided phase coherence is preserved over a sufficient large 
number of scattering events . In turn, these dis- 

coveries have cross-fertilized the field of light transport in 
the elastic regime [l3l fill ITEl fl(| . In this field, one of the 
hallmark of interference effects in elastic transport is the 
coherent backscattering (CBS) phenomenon |l7l ll8| : the 
average intensity multiply scattered off an optically thick 
sample is larger than the average background in a small 



angular range around the direction opposite to the ingo- 
ing light. This interference enhancement of the diffuse 
reflection off the sample is a manifestation of a two-wave 
interference. As such, it probes the coherence properties 
of the outgoing light and it has been extensively studied 
both experimentally and theoretically. It can be shown 
on general arguments, that the CBS enhancement factor 
(defined as the ratio of the backscattering CBS peak to 
diffuse background) never exceeds the value 2 and is ob- 
tained in the hclicity-preserving polarization channel for 
scatterers with spherical symmetry |19|. 

Whereas these interference modifications of transport 
are by now widely understood in the case of linear me- 
dia, recent experimental developments have required an 
extension of multiple scattering theory to the nonlinear 
case. Even if few studies already exists, they only cover 
the simpler case of classical linear scatterers embedded in 
a nonlinear medium pfl l2l| , whereas in our microscopic 
approach, the nonlinear behavior of randomly distributed 
scatterers will affect both the scattering processes and 
the average propagation. In particular, with the advent 
of laser cooling, on the one hand, it has become possible 
to study interference effects in multiple scattering of light 
by cold atoms |22t 12a. |2 ll 12a . l26j . In the regime where 
the saturation of the atomic transition sets in, atoms 
scatter light nonlinearly, i.e. the scattered light is no 
longer proportional to the incident one. One should note 
that important nonlinear effects are easily achieved with 
atoms even at moderate laser intensities. Considering 
a given driven optical dipole atomic transition, the or- 
der of magnitude of the required light intensity to induce 
nonlinear effects is given by the so-called saturation in- 
tensity /„ and is generally low. As typical examples, it 
is 1.6 mW/cm 2 for Rubidium atoms and 42 mW/cm 2 
for Strontium atoms, for their usual laser cooling tran- 
sitions. On the other hand, random lasers - mirrorless 
lasers where feedback is provided by multiple scattering 
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- have been realized experimentally [2^, |2^| . Here, 
nonlinear effects occur in the regime close to or above 
the laser threshold. Since, at least in the regime of co- 
herent feedback (2jj, interference is believed to play a 
decisive role in the physics of the random laser, a better 
understanding of the influence of nonlinearity (and am- 
plification) on the properties of coherent wave transport 
becomes necessary. 



II. MOTIVATION AND OUTLINE 

In a recent contribution |3lj , we have shown that non- 
linear scattering may fundamentally affect interference in 
multiple scattering. Indeed, in the perturbative regime of 
at most one scattering event with nonlinearity, there 
are now three (and no longer two) CBS interfering am- 
plitudes. Depending on the sign of the nonlinearity, i.e. 
depending whether nonlinear effects enhance or decrease 
the scattering cross section, the effect of this three-wave 
interference effect leads to a significant increase or de- 
crease of the nonlinear CBS enhancement factor. 

The purpose of the present paper is, on the one hand, 
to provide a detailed derivation of the equations for the 
nonlinear coherent backscattering signal used in 1311 . 
and, on the other one, to extend the treatment of [3l| 
to the case of atomic scatterers. Here, in contrast to 
the classical case, light is scattered inelastically, i.e. the 
scattered photons may change their frequencies. This 
leads to dephasing between interfering amplitudes and, 
consequently, to a reduction of the CBS enhancement 
factor in addition to the nonlinear modifications men- 
tioned above. Theoretical studies of this inelastic de- 
coherence mechanism have been so far restricted to the 
case of two atoms j33, |^ . Since the total (linear and 
nonlinear) elastic signal can be filtered out by means of a 
suitable frequency-selective detection, a clear experimen- 
tal study of inelastic, nonlinear CBS becomes possible. 
Please note that this would be otherwise very difficult to 
achieve since for weak intensities - the regime where our 
theory is valid - the linear signal generally largely domi- 
nates over the nonlinear one. In this paper, we will show 
that the enhancement factor for inelastically scattered 
light significantly exceeds the linear barrier two in cer- 
tain frequency windows. In contrast, the total enhance- 
ment factor - including also elastically scattered light - 
is diminished by nonlinear scattering. This is due to the 
negative sign of the total nonlinear component, since the 
total (elastic plus inelastic) scattering cross section is de- 
creased by saturation. 

The paper is organized as follows. In Sec. IHII we 
present the perturbative theory for nonlinear CBS of light 
scattered off a sample of cold two-level atoms. 'Perturba- 
tive' here means that we restrict ourselves to the regime 
of scalar - i.e., we forget the polarization of the photon 
- two-photon scattering with at most one nonlinear scat- 
tering event. This assumption is valid at sufficiently low 
probe intensities and not too large optical thicknesses. 



After shortly sketching the main results of the linear 
case, Sec. 1111 Al we derive equations for the nonlinear 
backscattering signal in Sec. IIII Bl The latter contains 
an inelastic and an elastic component. The latter again 
splits into a nonlinear and a linear part. In Sec. IIII CI 
supplemented by appendix [S] we show how to general- 
ize our scalar theory to the vectorial case by explicitly 
taking into account the light polarization degrees of free- 
dom. It is shown that nonlinear polarization effects lead 
to decoherence between interfering paths. In contrast 
to the linear case, this decoherence mechanism cannot 
be avoided by a suitable choice of the polarization de- 
tection channel. In order to emphasize the generality of 
our approach, we shortly discuss in Sec. IIII Dl a model of 
classical, nonlinear scatterers, which reproduces the elas- 
tic backscattering signal of the atomic model. In Sec. II VI 
we apply our theory to the case of a disordered atomic 
medium with slab geometry. We look at the dependence 
of the backscattering signal as a function of the optical 
thickness and of the detuning of the laser from the atomic 
resonance. In particular, we show that the enhancement 
factor for the inelastic component significantly exceeds 
the linear barrier two in certain frequency windows. Fi- 
nally, Sec. IVl concludes the paper. 



III. THEORY 

In this section, we present the perturbative theory for 
nonlinear coherent backscattering of light from a gas of 
cold two-level atoms. We first treat the linear component 
of the backscattering signal, which results from scatter- 
ing of independent photons. Thereby we introduce the 
reader, in Sec. IIII Al to standard methods used in lin- 
ear multiple scattering theory j3^, which we will then 
generalize to the non-linear case in Sec. IIII Bl 



A. Scalar linear regime 

1. One-photon scattering amplitude 

By definition, the linear component of the photo- 
detection signal is proportional to the incoming inten- 
sity, in particular to the number of photons in the initial 
laser mode. Since this implies that the photons are in- 
dependent from each other, it is sufficient to know how 
a single photon propagates in the atomic medium, see 
Fig. ^ This is equivalent to using the usual Maxwell's 
equations for a disordered medium |35j| . 

In the weak scattering regime, which we will consider 
throughout this paper, transport is depicted as a succes- 
sion of propagation in an average medium interrupted by 
scattering events. The important building block to prop- 
erly describe scattering and average propagation is the 
one-photon scattering amplitude by a single atom. For 
near-resonant scattering, and for atoms with no ground- 
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state internal Zeeman degeneracies, it reads : 

— Airi 



k(l - 2i8/T) ' 



(1) 



It can be derived from the elastically-bound electron 
model in the limit of small light detuning S = lu — oj at -C 
to, o; a t 5.5j. The atomic angular transition frequency is 
w a t whereas the atomic transition width T describes ra- 
diative decay. The photon wave number is k and the 
photon angular frequency is lu — ck (c being the vacuum 
speed of light). 

For simplicity we work here with scalar photons, i.e. 
we discard the vectorial nature of the light field. Scatter- 
ing is then fully isotropic and the differential scattering 
cross-section simply reads 
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where <jq is the on-resonance scattering cross-section. 

The scalar assumption is not a crucial one: as will 
be shown in Sec. IIII CI the following treatment can be 
generalized to the vectorial case. Please note however 
that the inclusion of internal degeneracies is not imme- 
diately simple and requires a separate treatment on its 
own. This is so because then the internal dynamics is no 
longer simple (optical pumping sets in). In this respect 
the results presented throughout this paper only apply to 
non-degenerate ground-state atoms. Please note also that 
internal degeneracies are already known to strongly re- 
duce the CBS effect in the linear regime 



2. Linear refraction index 



Between two successive scattering events occurring at 
r and r', the photon experiences an effective atomic 
medium with refractive index n^. Formally the resulting 
propagation is described by the average Green's function: 



G,(r,r') 



4tt r - 



where the refractive index is given by [36 



n u = 1 - 
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The imaginary part of n w describes depletion by scat- 
tering. This depletion gives rise to the exponential at- 
tenuation of the direct transmission through the sample 
(Beer-Lambert law) and defines, via the optical theorem, 
the linear mean-free path at frequency lu as 
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(6) 
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FIG. 1: Scattering path of a single photon entering the 
medium and leaving it in the backscattering direction to reach 
the detector. Straight lines depict average propagation in the 
effective medium while full circles depict scattering events la- 
beled by the r„. 



where AT denotes the density number of atoms in the 
sample. The weak scattering condition, where all the 
previous (and following) results are valid, then simply 
reads ki ^> 1. 



3. Linear radiative transfer equation 

We have now at hand all the necessary ingredients to 
write down the amplitude of a multiple scattering process 
as the one sketched in Fig. ^ We consider a scattering 
volume V exposed to an initial monochromatic field with 
amplitude Eq propagating along axis z. The transverse 
area of the scattering volume is S. Since ki ^> 1, a 
semi-classical picture using well-defined scattering paths 
is appropriate. For a given scattering path C„ = (ri — > 
• • • — > r n ) labeled by the collection of scattering events, 
the corresponding far-field amplitude radiated at position 
R of the detector placed in the backscattering direction 



S(C n ) 



AnkR 



A(C n )E . 



(7) 



The complex amplitude A(C n ) is simply a product 
of one-photon scattering amplitudes 0} and of Green's 
functions (@J: 



'n-X 



A(C n ) = kS^e lkn ^ +z ^ 11 S* w G w (ri,r i+1 ) (8) 



where is the distance at which scattering event i occurs 
from the boundary of the medium. The superposition 
principle then gives the total electric field amplitude £ 
as a sum over all possible scattering paths C n : 



£ = - 
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A = Y,A(C n ) (9) 
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The total average intensity is obtained by squaring 
and averaging over all possible scattering events. We 
define the total dimensionless bistatic coefficient as: 

7^ - Hjr <|£| 2 )di S .av. = (\A\ 2 ) dis . av . (10) 

We now assume complete cancellation of interference 
effects between different scattering paths (random-phase 
or Boltzmann approximation) . We then obtain the back- 
ground (or 'ladder') component of the backscattering sig- 
nal: 

7^ * = £ ' / dv x . . . dr n \A(C n )\ 2 (11) 

n=l JV 

This formula has a well-defined limit when E — > oo and 
thus can be applied to slab geometries. Please note that, 
in writing 1|11[1 , we have also discarded recurrent scatter- 
ing paths, i.e. paths visiting a given scatterer more than 
once. Both approximations are justified in the case of a 
dilute medium, hi ^> 1 37] . 

We rewrite (|1H as 

^=/|^(r)e-^ (12) 

with 

oo „ 

I u (r) = tT'l 1 + V Af n / dv 1 . . . dr n 

n=i Jv 

n 

e^ l/£ ni^ G -( r i' r i+i)| 2 ' ( 13 ) 

1=1 

where r n+ i = r. This dimensionless function describes 
the average light intensity at r, in units of the incident 
intensity I = cqcEq/2 (in W/m 2 ) with eo the vacuum 
permittivity. The first term in Eq. (|13[) represents the 
exponential attenuation of the incident light mode, i.e. 
light which has penetrated to position r without being 
scattered (Beer-Lambert law). The remaining term de- 
scribes the diffuse intensity, i.e. light which has been 
scattered at least once before reaching r. From Eq. (|13fl . 
one can easily show that 1^ (r) fulfills the radiative trans- 
fer integral equation Q: 

J u (r)=e-*/' + y ^dr'|G w (r,r')| 2 ^(r') (14) 

The required solution of Eq. I|14[l can be obtained numer- 
ically by iteration starting from I w (r) = 0. 

4- Linear CBS cone 

In fact, the preceding Boltzmann approximation 7^ « 

L^V is wrong around the backscattering direction. In- 
deed, on top of the background ladder component, one 



observes a narrow cone of height C el and angular width 
In the regime hi ^> 1, this so-called 
CBS cone arises from the interference between ampli- 
tudes associated to reversed scattering paths C„ = (ri — > 
■ ■ • — > r n ) and C„ = (r n n). Of course single 

scattering paths where n = 1 do not participate to this 
two- wave interference (since they are exactly identical to 
their reversed counterparts) and must be excluded from 

. Thereby, we obtain the interference (or 'crossed') 
contribution as: 

C ci 1, = E^ J v dv l ...dv n A(C n )A*{C n ) (15) 

Thus, the bistatic coefficient in the backscattering di- 
rection reads 7^ = L^' + cj^' . From Eq. ©, we verify 

that the reciprocity symmetry A(C n ) = A(C n ) is fulfilled 
for scatterers without any internal ground-state degen- 
eracies. This allows us to rewrite Ijl5|l as 

<*?-/§ (Mr) -.-*«)•-"' 

where is the single scattering contribution. Hence, 
the linear CBS enhancement factor, defined as 

^=1 + ^/^=2-3$/^, (17) 

is always smaller than two. It equals two if single scat- 
tering can be filtered out, see Sec. IIII Cl 

B. Scalar nonlinear regime 

At higher incident intensities, the successive photon 
scattering events become correlated. Indeed absorption 
of one single photon brings the atom in its excited state 
where it rests for a quite long time T _1 without being able 
to scatter other incident photons. This means that sat- 
uration of the optical atomic transition sets in, inducing 
nonlinear effects and inelastic scattering. In a perturba- 
tive expansion of the photo-detection signal in powers of 
the incident intensity, the leading nonlinear term arises 
from scattering of two photons. In order to generalize the 
above linear treatment to the two-photon case, we first 
need to remind some relevant facts about scattering of 
two photons by a single atom [32I ]. 

1. One-atom two-photon inelastic spectrum 

The two-photon scattering matrix S contains an elastic 
and an inelastic part. The elastic part corresponds to two 
single photons scattered independently from each other, 
whereas the inelastic part describes a 'true' two-photon 
scattering process, where the photons become correlated 
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and exchange energy with each other. To obtain the in- 
tensity of the photo-detection signal, the electric field 
operator E (evaluated at the position of the detector) is 
applied on the final two-photon state |/) = S\i), with 
\i) the initial state. Since E annihilates one photon, this 
yields a single-photon state = E\f), which describes 
the final state of the undetected photon. Like the scat- 
tering matrix S, it consists of an elastic and an inelastic 
component: 



(18) 



The inelastic part |"0in) is a spherical wave emitted by the 
atom, whereas the elastic part |^ e l) is a superposition of 
scattered and unscattered light, thereby taking into ac- 
count forward scattering of the undetected photon. (For- 
ward scattering of the detected photon does not need to 
be taken into account, since the detector is placed in the 
backscattering direction.) Finally, the norm / = (?/>|?/>) 
of \ip) defines the intensity of the photo-detection signal. 
According to Eq. (|18|l . / is the sum of the following three 
terms: 



= (^cll^el), 
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(19) 

(20) 
(21) 



So far, everything is valid for any two-photon scattering 
process with an elastic and an inelastic component. In 
the specific case of a single atom, the following result is 
obtained: 
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with the incident intensity Iq, and the saturation param- 
eter s defined by [38| : 
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where d is the atomic dipole strength and I s the satura- 
tion intensity of the atomic transition. 

The first term, Eqs. I|19l22fl . which arises from two 
photons scattered independently from each other, repro- 
duces the linear single-photon cross section Air a = |5 W | 2 , 
see Eq. ([Q). The following two terms correspond to non- 
linear elastic and inelastic scattering, respectively. For 
the case of a single atom, the perturbative two-photon 
treatment is valid for s <§; 1, i.e. if the nonlinear terms 
are small compared to the linear one. 

The frequency spectrum of the elastically scattered 



light is simply F c \{uj') = (ijp + i^) <5(o/ — uj) whereas 
the frequency spectrum of the inelastically scattered light 
is Fi n (w') = 7; n -P(w'). The continuous spectrum P(uj') 



r(2)> 




A 



FIG. 2: In the perturbative approach, we assume a single non- 
linear two-photon scattering event (□), but arbitrarily many 
linear scattering events (•). One of the two photons is fi- 
nally annihilated by the detector, thereby defining the photo- 
detection signal, whereas the other one is scattered into an 
arbitrary direction. 



is normalized to unity according to J du}'P(u)') = 1. It is 
obtained as follows 321: 



1 



1 



iT/2 25 - 5' + iT/2 



(26) 



where 8' = u>' — oj at denotes the final detuning. This 
inelastic spectrum consists of two peaks with width T, 
one located at the atomic resonance (u/ = w a t), and 
the other one twice as far detuned as the incident laser 
(a/ = w at + 25). For 6 < T/2, the two peaks merge to 
a single one centered at u>' = u>. Please note that, by 
going beyond the two-photons scattering approximation, 
one would then get three peaks as predicted by the non- 
perturbative calculation of the inelastic spectrum, also 
known as the Mollow triplet [38| . 

2. Nonlinear scattering in a dilute medium of atoms 

Now, we generalize the above single-atom treatment to 
a multiple scattering process in a dilute medium of atoms. 
First, we note that the above perturbative treatment - in 
particular Eqs. (|19I21|) - remains valid for any form of 
the scattering sample, be it a single atom, two atoms, or 
arbitrarily many of them. An important difference from 
the single-atom case, however, is that the total weight 
of nonlinear processes may be drastically enhanced if the 
sample has a large optical thickness b = L/£, where L 
is the typical medium size. This implies that the condi- 
tion s <C 1 is not sufficient to guarantee the validity of 
the perturbative approach. Instead, as we will argue in 
Sec. IIVI the perturbative condition reads sb 2 <C 1. 

A typical two-photon scattering path is sketched in 
Fig- HI Here, the incoming photons propagate at first in- 
dependently from each other to position r inside the dis- 
ordered atomic medium, where they undergo a nonlinear 
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scattering event. One of the two outgoing photons then 
propagates back to the detector. The possibility that the 
two photons meet again at another atom can be neglected 
in the case of a dilute medium, similar to recurrent scat- 
tering in the linear case [3?]]. We can hence restrict our 
analysis to processes like the one shown in Fig. [21 with 
arbitrary numbers of linear scattering events before and 
after the nonlinear one. Thus one of the two incoming 
photons undergoes n > elastic scattering events (la- 
beled by Uj), while the other undergoes m > elastic 
scattering events (labeled by vj), before merging at r 
where they undergo the inelastic scattering event. One 
of the outgoing inelastic photons reaches back the detec- 
tor after having undergone I > elastic scattering events 
(labeled by positions Wk). For the other undetected in- 
elastic photon, we may assume, without any loss of gen- 
erality, that it does not interact anymore with the atomic 
medium. This interaction would be anyway described by 
a unitary operator (as a consequence of energy conserva- 
tion), which does not change the norm of the state of 
the undetected photon defining the detection signal. 

In general, the state of the inelastic undetected pho- 
ton corresponding to a scattering path C defined by the 
position r of the two-photon scattering event and by the 
collection of positions of all one-photon scattering events 
C = {u, v, r, w} is given as follows: 

|-0 in (C)) = e «»"t*«i+^i) 

n m 

x JJS , w G w (u i ,u i+ i)]^[5 w G w (Vj,Vj + i) 



i=i 



J U 1 



k=l 



e tfcn «' z »i x < (27) 
12 n > or m > 
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FIG. 3: The elastic component |^>(r, {u, v, w}) e i) of the un- 
detected photon's state arises from a superposition of the fol- 
lowing three processes: (a) elastic scattering of both photons 
at r, (b) only the detected photon scattered at r, and (c) only 
the undetected photon scattered at r. The last two diagrams 
are necessary to take into account the nonlinear average prop- 
agation of the undetected (b) or detected (c) photon. 



contained in the calculation of the inelastic component. 
In contrast to the single-atom case, however, the elastic 
component will enter in the calculation of the nonlinear 
average propagation, i.e. the nonlinear modification of 
the refractive index (Kerr effect), and will be discussed 
later. At first, we concentrate on the processes of non- 
linear scattering, i.e. processes changing the direction of 
propagation of the detected photon. 

As for the linear case, we still assume the same di- 
lute medium approximations to hold for the 'ladder' and 
'crossed' contributions. Thus, in order to calculate the 
average photo-detection signal, we just keep scattering 
diagrams obtained by reversing the path of the detected 
photon. Furthermore, we also neglect interference be- 
tween diagrams where the nonlinear scattering event oc- 
curs at different atoms. This is justified in the dilute case 
since the overlap between two spherical waves emitted at 
r and r' vanishes if fc|r — r'| ^> 1. 



3. Nonlinear ladder contribution 



with u n+ i = v m+ i = wi + i = r, IT/ the projector on 
photons states at frequency w' and |^>i n ) the inelastic final 
state of the one-atom case, Eq. (JTSJ. Since the inelastic 
two-photon scattering event takes place at position r, this 
state describes an outgoing spherical wave emitted at r. 
Furthermore, note that if the two incoming photons do 
not originate both from the incident mode, i.e. if n > or 
m > 0, a factor 2 arises due to the fact that the incoming 
photons can be distributed in two different ways among 
the paths {u} and {v}. 

The elastic component \ip c i(C)) is obtained in a sim- 
ilar way. However, as in the single-atom case, we must 
take into account forward scattering of the undetected 
photon, at the position r of the nonlinear event. This is 
done by considering the superposition of two diagrams 
where the undetected photon is scattered or not scat- 
tered at r, see Fig. |3fa,b). Since this approach exactly 
parallels the one known from the single-atom case , it 
is unnecessary to present the complete calculation of the 
elastic component in detail - all relevant ingredients to 
perform the generalization to the multi-atom case will be 



To obtain the inelastic component of the average 
backscattering signal, we first get the total final state 
of the undetected photon by summing Eq. I|27|l over all 
possible different scattering paths C. Then we insert 
this result into Eq. I|21|l and we finally average over the 
random positions of the scatterers. As argued above, 
only identical or reversed scattering paths are retained 
in the average, giving rise to the background ('ladder') 
and interference ('crossed') component. Thus, the inelas- 
tic background component reads as follows: 

r v xfn+m+l + l r n m 1 

* e W/ V n*n^n^ 

JV (n,m,0=0 JV i=l 3=1 fe=l 

GMC)|v*n(c)> x \\ if " = ™ = o (28) 

I 1/2 otherwise 

Note that some care must be taken not to sum twice 
over the same scattering path. In particular, any ex- 
change of the two incoming parts {u} and {v} leaves the 



total scattering path unchanged since the two incoming 
photons are identical. For this reason, a factor 1 /2 must 
be inserted at the end of Eq. I)28[l. Again, as in Eq. (|27|1 . 
the case n = m — is exceptional, since then there is 
no elastic scattering events before the nonlinear one: the 
two incident photons remain in the same mode. 

If we insert now Eq. I|27l) into Eq. (|28|1 . we simply ob- 
tain the inelastic nonlinear ladder contribution as: 



-(2) 



dr 



(2/2(r)-e- 2 ^) f dJP{w')lAr) 



(29) 

with Iu(r) the linear average intensity, see Eq. Q14JI. In 
order to interpret this result, we first note that the in- 
elastic intensity radiated by the atom at position r is 
proportional to the mean squared intensity at r. An al- 
ternative, physically transparent derivation of the latter 
can be performed as follows: we write the local field am- 
plitude A — exp(— z/2£) + Ad as a sum of coherent and 
diffuse light amplitudes. The latter term exhibits a Gaus- 
sian speckle statistics [3!j, i-e. (Re„4o) = (Invto) = 0, 
2((ReA D ) 2 ) = 2((ImA D ) 2 ) = (\A D )\ 2 ) and (|^| 4 ) = 

2(|-4d| 2 ) 2 - Thereby, we obtain for the mean squared 
intensity: 



(\A\ 4 ) = e- 2 ^ + (\A D \ 4 )+4e-^(\A D \ 2 ) (30) 
= 2(|^| 2 > 2 -e- 2 ^. (31) 



Inserting the average intensity, I u — (l^l 2 ), we immedi- 
ately recognize the first integrand in Eq. (|29() . Then, the 
atom emits a photon with frequency distribution P(u>'). 
Finally, due to time reversal symmetry, the propagation 
of this photon from r to the detector is described by the 
same function I u i (r) which represents propagation of in- 
coming photons to r. 

Concerning the elastic component, the diagrammatic 
calculation via Eq. l(2"U)l. see also Fig.[3{a,b), shows that 
the above argument can be repeated in the same way 
- except for the fact that the detected photon does not 
change its frequency. Furthermore, a factor —2 is taken 
over from the single-atom expression, cf. Eqs. (|23I24|I . 
Thereby, we obtain: 



L 



(2,scatt) 



^/g(2/ 2 (r)-e- 2 ^)/.(r). (32) 



The index 'scatt' reminds us that we have treated only 
nonlinear scattering so far. Below fSec. 1111 B 5|) . we will 
add nonlinear average propagation, which contributes to 
the clastic nonlinear component, too. 



4- Nonlinear crossed contribution 






^37 



VZ7 



FIG. 4: In the presence of nonlinear scattering (□), there 
may be either (b) two, or (c) three interfering amplitudes con- 
tributing to enhanced backscattering, apart from single scat- 
tering (a), which only contributes to the background. In gen- 
eral, the case (c), which corresponds to maximum enhance- 
ment factor three, is realized if either both incoming photons, 
or one incoming and the outgoing detected photon exhibit at 
least one linear scattering event (•) besides the nonlinear one. 



admits more possibilities to reverse the path of the de- 
tected photon. This is due to the photon exchange sym- 
metry at the nonlinear scattering event, which does not 
allow to distinguish which one of the two incoming pho- 
tons finally corresponds to the detected or undetected 
one. As evident from Fig. 0Jc), each multiple scattering 
path where both incoming photons, or one incoming and 
the outgoing detected photon, exhibit at least one lin- 
ear scattering event besides the nonlinear one, has two 
different reversed counterparts, leading in total to three 
interfering amplitudes. 

If we look at the scattering process shown in Fig. El the 
two reversed counterparts are obtained by exchanging the 
outgoing detected photon {w} with either one of the in- 
coming photons {u} or {v}. Since both cases are identical 
in the ensemble average, we may restrict ourselves to one 
of them, let us say {v}. We thus note C = {u, w, r,v} 
the reverse path corresponding to C = {u, v, r, w} when 
{v} and {w} are exchanged. In total, we obtain for the 
inelastic interference component: 



(2) 



dr £ 



TV 



■n+m+l+l 



v 47Tfc 2 S 

(n,m,l)= 



n m 



n dm ii n dw k (Vin(c)iv>i„(c)> 

if m = I = 



V i=i j=i k=l 



1 1 otherwise 



It remains to calculate the 'crossed' contribution, i.e. 
interference between reversed paths. In contrast to the 
linear case, where there are always two interfering ampli- 
tudes (apart from single scattering), the nonlinear case 



Here, the case m = I — identifies processes where the 
two reversed paths C and C are indistinguishable. Set- 
ting their contribution equal to zero accounts in partic- 
ular for the single scattering case depicted in Fig. EJa), 
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i.e. n = m = I = 0, which does not contribute to the 
interference cone. The case Fig. E{b) remains with two 
contributions [n = m = 0, I > 0, and n = I = 0, m > 0, 
respectively) in Eq. corresponding to the fact that 

two amplitudes interfere. Finally, the case (c) of three 
interfering amplitudes is reflected in Eq. ()33|l by the ab- 
sence of the exchange factor 1/2, as compared to the 
background, Eq. 1)28)1 . Thereby, the interference contri- 
bution can, in principle, become up to two times larger 
than the background. 

If we insert the state of the undetected photon, 
Eq. (|27|) . into Eq. l|33[l . we encounter the following ex- 
pression 

OO p 

g u M T ) = e ^ n * + Vat / «&•!...£&•» 

n=l J V 

n 

e*^^*l[S u G u (r u r t+1 )S^GZ,(r u r t+1 ), (34) 
i=i 

which generalizes the local intensity, Eq. I|13|l , to the case 
where two different frequencies occur in the interfering 
paths. Numerically, it can be obtained as the iterative 
solution of: 



— e ik(r 



Jv 



(35) 



This function describes the ensemble-averaged product 
of two probability amplitudes, one representing an in- 
coming photon with frequency uj propagating to position 
r, and the other one the complex conjugate of a photon 
with frequency uj' propagating from r to the detector. If 
uj 7^ uj' , then these amplitudes display a nonvanishing 
phase difference both due to scattering and to average 
propagation in the medium. This leads on average to a 
decoherence mechanism and consequently to a loss of in- 
terference contrast. Indeed, both the complex scattering 
amplitude, Eq. (JIJ, and the refractive index, Eq. ©, de- 
pend on frequency. In contrast, the phase difference due 
to free propagation (i.e. in the vacuum) can be neglected 
if T£ <gC c . w hich is fulfilled for typical experimental pa- 
rameters [4(J, |41j . In the case uj = uj' of identical frequen- 
cies, (? w , w (r) = / w ( r ) reduces to the average intensity, see 
Eq. (nj). 

In terms of the iterative solution of Eq. 134|) . the in- 
elastic interference term, Eq. I|33ll . is rewritten as follows: 



=4, / du>'P(v') I ^ / w (r)| 5w , w ,(r)| 



f/r 



-e-^Re{ e ^-»:') fo ^ iW ,(r)} 



(36) 



with £' the linear mean free path at frequency uj' . In the 
elastic case, where u>' = u>, dephasing between reversed 





FIG. 5: Diagrammatic description of nonlinear propagation. 
A two-photon process (solid lines) interferes with two inde- 
pendent single photons (dashed lines). Only one of the latter 
(the undetected photon) is scattered at □, thereby modify- 
ing the propagation of the detected photon (a) between two 
scattering events at positions 1 and 2, (b) on the way to the 
detector after the last scattering event at position 1, or (c) in 
the coherent mode before the first scattering event at position 
1. 



scattering paths does not occur, and the expression ll3f)|) 
simplifies to: 



C, 



(2,scatt) 
J 



-8s 



dr 



(/ w (r) 3 - 2I u (r)e 



■2z/e 



+ e 



-3z/e 



(37) 

Since, in the elastic case, there is no loss of coherence 
due to change of frequency, the elastic interference com- 
ponent, Eq. 13 7|) . is completely determined by the relative 
weights of the one-, two-, and three- amplitudes cases ex- 
emplified in Fig. This can be checked by rewriting the 
background and interference components, Eqs. I|32I37|I . 
in terms of diffuse and coherent light, respectively, i.e. 
by writing / = Ijj + cxp(— z/t). One obtains: 

L (2,scatt) k ^ 3z/e + Moe -2 Z /i + Q^-z/i + 2/ 3 ^ 



(2.scatt) 
el 



OC 



41 D e 



-2zl 



\2I 2 D e 



(a) 



-z/l 



(c) 



(38) 

where the brackets denote the integral over the volume V 
of the medium, and (a,b,c) correspond to the three cases 
shown in Fig. 0] identified by different powers of diffuse 
or coherent light. As expected, the three-amplitudes case 
(c) implies an interference term twice as large as the back- 
ground. In the two-amplitudes case (b), a small compli- 
cation arises, since one of the two interfering amplitudes 
is twice as large as the other one (i.e. the one where both 
incoming photons originate from the coherent mode), cf. 
the discussion after Eq. 127|l . In this case, the interfer- 
ence contribution, 2x1 + 1x2 = 4, is smaller than the 
background, 2x2 + 1x1 = 5. Finally, as it should be, 
the single scattering term (a) is absent in the interference 
term C, 



(2,scatt) 
ol 



5. Nonlinear average propagation 

So far, we have only considered processes of nonlin- 
ear scattering where the direction of propagation of the 
detected photon is changed. It remains to take into ac- 
count nonlinear average propagation. This is described 
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by those processes where, in one of the two interfering 
amplitudes, the detected photon is not scattered at the 
position r of the nonlinear event |44j . The corresponding 
diagrams are depicted in Fig. [5] where the two interfering 
amplitudes are represented by the solid and dashed lines, 
respectively. Here, the solid lines correspond to an inelas- 
tic two-photon scattering process (like the one shown in 
Fig- 13 j whereas the dashed lines represent an elastic pro- 
cess, where the two photons are independent from each 
other, see Fig-Ofc). Hence, their interference contributes 
to the nonlinear elastic component of the photo-detection 
signal, cf. Eq. (|2Ull . 

The three diagrams shown in Fig. [5] differ only by the 
fact that the nonlinear propagation event takes place ei- 
ther between two scattering events at position 1 and 2 
(a), on the way to the detector, i.e. after the last scat- 
tering event at position 1 (b), or in the coherent mode, 
i.e. before the first scattering event at position 1 (c). At 
first, let us examine the case (a). We imagine that each 
of the three dots • may represent an arbitrary number of 
scattering events. [Only note that the number of events 
corresponding to the dots and 1 and 2 must be larger 
than zero - otherwise, the diagram Fig. EJa) would be 
identical to Fig.^b) or (c).] According to the theory of 

HTTA1 



linear radiative transfer outlined in Sec. lIII Al the ladder 
diagrams corresponding to the two incoming photons ar- 
riving at 1 (position ri) and at the nonlinear event □ 
(position rs) yield the linear local intensities /^(ri) and 
I u (ra), respectively Likewise (due to reciprocity sym- 
metry) , the propagation of the outgoing detected photon 
from 2 (position r 2 ) to the detector - with arbitrary num- 
ber of scattering events in between - is given by I ul (r2). 
Hence, the only ingredient which we have to calculate is 
the nonlinear propagation between 1 and 2. Note that, 
when taking the average over the position r3 of the non- 
linear event, non- negligible contributions arise only if r3 
is situated on the straight line between ri and v^, since 
this is the only way to fulfill a stationary phase (or phase 
matching) condition. Thereby, the 'pump intensity' en- 
tering in the nonlinear propagation is given by the av- 
erage value of the local intensity on this line, which we 
denote by (/ w )n^r 2 - We do not want to present the 
complete calculation here (this requires to calculate at 
first the case of a single atom, which can be done with 
the techniques described in 32]), but just give the final 
result: 



G^\v^v 2 ) =|G w ( ri ,r 2 )| ! 



2sris 



(39) 



From this, we deduce the following value for the nonlinear 
mean free path: 



(40) 



which is consistent with Eq. I|39[) . if we expand the re- 
sulting propagator (where the mean free path appears in 
the exponent) up to first order in s. The same result is 
also obtained in the case of diagram Fig. EJb), i.e. for 



the propagation after the last scattering event. Hence, 
the corresponding propagator (first order in s) reads: 



GL nl < b) (ri) 



, 2sz\ 



(I 



(41) 



where ro = 1*1 — zie z , with e z the unit vector pointing 
in the direction of the incident laser, denotes the point 
where the photon leaves the medium. In the case (c), 
a small complication arises since the photons arriving at 
the nonlinear event □ may originate both from the coher- 
ent mode, which reduces the two-photon scattering am- 
plitude by a factor 1/2, cf. the discussion after Eq. I|27|l . 
Hence, the nonlinear mean free path for photons from 
the coherent mode reads: 



1 



j(l-2sUr) + S e- z l 



with the corresponding propagator 



GL nl < c) (n) 



= e-^ffl (2L(r) _ ( 



t /t\ 

I r 



(42) 



(43) 



The difference between the mean free paths, Eqs. (|40I42II . 
can also be understood as a consequence of the different 
properties of intensity fluctuations for diffuse and coher- 
ent light, see Eq. (|30|) - which determine the nonlinear 
atomic response. 

In total, we obtain for the background component: 



-(2, prop) _ 



el 



4,(ri) 



G, ( ?'' b) ( 



(44) 



In the case of a slab of length L, Eq. I|44|) can be simplified 
to: 



£(2,prop) _ 9 



—I u (z) (2Uz) 2 



-2z/t _ e -z/A 



(45) 



Concerning the interference component, we find the 
same phenomenon which we have already observed in the 
case of nonlinear scattering: if we exchange outgoing and 
incoming propagators, we find twice as many 'crossed' as 
'ladder' diagrams, see Fig- El I n particular, the diagrams 
(d,e,f), which could be seen as a modification of the lin- 
ear refractive index by the local crossed intensity - thus 
affecting the (ladder) average propagation - are not con- 
sidered in previously published papers, concernin g ei ther 
classical linear scatterers in a non-linear medium |2(1 l2ll ] 
or nonlinear scatterers in the vacuum |3lt |42j . Even if, 
at first sight, these diagrams look unusual, our numerical 
calculations (see section HlID|) suggest that they play an 
important role, at least in our situation where nonlinear 
scattering and nonlinear propagation originate from the 
same microscopic process. 
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a) 



A 





nents discussed above: 



C = C, 



-(i) 
(i) 



- (2,scatt) 
(2,scatt) 



r (2, prop) r (2) 
L el +^in > 



c, 



c. 



(2, prop) 
cl 



c 



(2) 



(48) 
(49) 



Before we present the numerical results in Sec. II VI we will 
generalize the above results to the vectorial case. This is 
important since polarization does not only lead to slight 
modifications for low scattering orders, as in the linear 
case. Apart from that, we will see that it also induces de- 
coherence between reversed paths, thereby reducing the 
nonlinear interference components. 



FIG. 6: Interference contributions from nonlinear propaga- 
tion. The diagrams (a-f) are obtained from the ladder di- 
agrams, see Fig. |5J by reversing the paths of the respective 
photons. Just as in the case of nonlinear scattering, there are 
twice as many diagrams contributing to the interference cone 
than to the background. 



C. Incorporation of polarization : vectorial case 

First, including the polarization modifies the scalar ex- 
pressions, Eqs. 11 Kill , for the linear mean free path and 
the atom-photon scattering amplitude by a factor 2/3: 



Due to the reciprocity symmetry (remember that non- 
linear propagation contributes to the elastic component, 
i.e. no decoherence due to change of frequency ), each of 
the diagrams in Fig. Ogives the same contribution as the 
corresponding ladder diagram in Fig. [SJ Hence, to first 
approximation, the interference contribution from non- 
linear propagation equals twice the background, Eq. (|44Jl . 
Some care must be taken, however, if photons arriving at 
(or departing from) the nonlinear event (or position 1) 
originate from the coherent mode. In such cases, it may 
happen that some of the diagrams depicted in Figs. Eland 
El coincide, and we should not count them twice. [This 
is analogous to the distinction between the cases (a,b,c) 
in Fig. 0] or to the suppression of single scattering in the 
linear case.] 

Taking this into account [for details, we refer to the 
discussion after Eq. (|A8|) in appendix |A"]. we find: 



£f(2,prop) _ 2^(2, prop) 



+I(ri)e" 



dr 1 

r Yi 

-*i/e Q 



G (nl,c) ( 



(46) 



In the case of a slab, we obtain: 



C. 



(2, prop) 
ol 



2L 



(2, prop) 
il 



-z/t -25 

e ' — e 



: L dz , s 
:]s ' T I{z) 



(47) 



where b = L/t denotes the (linear) optical thickness of 
the slab. 

Thereby, we have completed the perturbative calcula- 
tion of the backscattering signal for the scalar case. The 
total signal is obtained as the sum of the various compo- 



S/jj — 



— 67ri 



k 2 



fc(l - 2id/T) 



(50) 
(51) 



The Green's function, Eq. (J3J, remains unchanged, ex- 
cept for the fact that the modified expression for the 
mean free path, Eq. I|50|) . must be inserted in the re- 
fractive index. However, the angular anisotropic charac- 
ter of the atom-photon scattering is not yet contained in 
Eq. 1)5 l[l. This is treated by projection of the polarization 
vector as follows: if the photon, with incoming polariza- 
tion ei, is scattered at rj, and the next scattering event 
takes place at r-2, the new incoming polarization reads: 



£ 2 — A ri>ra ei, 



(52) 



where A ri!l . 2 denotes the projection onto the plane per- 
pendicular to ri — t-2 . Finally, the detection signal after n 
scattering events is obtained as e* D ■ e„, with the detector 
polarization e^. 

Thus, the linear background ('ladder') contribution 
reads, cf. Eqs. ipTTJ) : 



— ^A/""" 1 / dv 2 . 
Al Jv 



.dr n e 



\\ |^Gc(ri,r i+1 )| 2 



-Zn/t 



■ A ri ,r 2 e L 



(53) 



where denotes the initial laser polarization. By choos- 
ing a given circular polarization, for example = 
0)/v2, and by detecting the signal in the helicity- 
preserving h\\h polarization channel {eu = e* L ), then the 
single scattering contribution in Eq. I|53|l (n — 1 term) is 
filtered out. We thus recover the enhancement factor 2, 
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a) 



FIG. 7: Polarization vectors associated to the two-photon 
scattering matrix for two reversed scattering amplitudes (a) 
and (b). Note that the corresponding reversed scattering am- 
plitudes, Eqs. 1541561 are different - even in the helicity pre- 
serving polarization channel, i.e. if £2,3 = e|,3- This leads to 
a reduction of the CBS interference cone by a factor 3/4, on 
average. 






FIG. 8: Polarization dependence of nonlinear propagation for 
(a) ladder and (b,c) crossed diagrams. 



meaning = . Apart from that, however, polar- 
ization does not play a very important role: the distri- 
bution of higher scattering orders n > 1 is only slightly 
modified, and the reciprocity symmetry remains valid, 
provided that en = e* L . 

The situation changes in the nonlinear regime of two- 
photon scattering. With the initial and final polariza- 
tions ei t 2 and £3,4, respectively, see Fig. 0a), the po- 
larization dependent term of the two-photon scattering 
matrix reads: 

Sp=l[(ei-el)(e 2 -e* 3 ) + (e 1 -e* 3 )(e 2 -et)]. (54) 

The prefactor 1/2 is chosen such that S p represents cor- 
rectly the polarized scattering amplitude in units of the 
corresponding scalar one. From Eq. Q54JI. the photon ex- 
change symmetry becomes evident: the outgoing photon 
3, e.g., can equally well be associated with the incoming 
photon 1 or 2. If we trace over the undetected photon, 
which we may label as photon 4, for example, we obtain 
for the ladder component: 



nW( ei ,e 2 ; e3 ) = El 5 pl 



l £ 2 • e 3 l 



I * |2 

ki • e 3 l 



+2Re{( ei -e|)(e 2 -e^)(e 3 -et)}] . (55) 



If we assume a random uniform distribution for the polar- 
ization vectors, we obtain (Il( L )) = 2/9, which is smaller 
than the linear counterpart (|e„ • e* D \ 2 ) = 1/3. Hence, 
in the vectorial case, the relative weight of the nonlinear 
contribution is approximately one third smaller than in 
the scalar case - at least far inside the medium, where 
the polarization is sufficiently randomized. 

Concerning the interference ('crossed') contribution, 
we exchange the direction of the outgoing detected pho- 
ton 3 and of one of the incoming photons, for example 
photon 2. Note that we obtain in general different po- 
larizations £2,3 for the reversed counterparts of £2,3, see 
Fig. |7fb) . Indeed, the reversed photons have the same 
polarizations, £2,3 = £2,3, only if the laser and detector 
polarizations are identical (f£> = £l). Consequently, the 



scattering amplitude for the complex conjugate photon 
pair reads: 



S P = \ [(ei ' e* 4 )(e 3 • e* 2 ) + I 



£1 ' £2) • el)] 



(56) 



Note that even in the helicity-preserving polarization 
channel, i.e. £2.3 = £23, the reversed scattering ampli- 
tudes, Eqs. H54I56(I . are in general not equal. Only the 
first term, where photon 2 is associated with photon 3, 
remains unchanged if those two photons are reversed. As 
a consequence, the polarization induces a loss of coher- 
ence, i.e. a reduction of the crossed term as compared to 
the scalar case. The sum over the polarization of photon 
4 yields: 



n^( ei . 



£2, £3, £3, £2) 



= i [(£2-£5)(£2-£l) + (£2-£^)(£l-£l)(£-2-£D 
+ ■ £^)(£ 2 • £l)(£2 • el) + ( £l • £*)(£ 2 ■ ■ £~3*)] . 

(57) 

If we assume £2,3 = £23, i-e. the h \\ h channel, we 
obtain ttl( c )) = 3/18 on average. Hence, in this channel, 
the polarization-induced loss of contrast is approximately 

(n( c V(n (i) ) = 3/4. 

Finally, to obtain the polarization dependence of non- 
linear propagation, we label the photons as shown in 
Fig- HI Let us first examine the ladder term, Fig. [Sfa). 
The solid lines are described by the two-photon ampli- 
tude, Eq. I]54p. whereas the dashed lines give the complex 
conjugate of (£2 • £4) (ei • £3) . After integration over photon 
4, the result is 

n( i 'P I '°P)(£ 1 ,£ 2 ,£ 3 ) =~ [(£1 • e* 2 )(e 2 ■ e*)(e 3 • ej) 

+ (£1 • £s)(£2 • e* 2 )(e 3 ■ c?)] . 

Concerning the crossed diagrams, we distinguish between 
the two cases shown in Fig.[H{b) and (c). [In Fig. El this 
corresponds to (a,b,c), on the one hand, and (d,e,f), on 
the other hand.] As for the case (b), nothing changes 
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since the reversed photon does not participate in the non- 
linear event. In case (c), we obtain: 

n( c >P r °P)( ei ,e 2 , e 3,e2,e3) = \ [(« • e|)(e~ 3 ' ^*)(e 3 • ej) 

+ (e 1 -e 2 *)(e~ 3 -e|)(e 3 -eD]. (59) 

When determining the average values of the nonlinear 
propagation terms, it must be taken into account that ei 
and 63 are not independent from each other, since they 
propagate in the same (or opposite) direction. Thus, we 
find (n( L 'P r °P)) = 1/3 and (n( c >P ro P)) = 1/6. Hence, the 
loss of contrast equals 1/2 in case (c), whereas reciprocity 
remains conserved (i.e. no loss of contrast) in case (b). 
Averaging over (b) and (c), this yields the same contrast 
3/4 as for nonlinear scattering. 

What remains to be done to obtain the vectorial 
backscattering signal is to incorporate the above expres- 
sions into the corresponding scalar equations. The re- 
sulting equations can be found in appendix lAl together 
with a description of the Monte-Carlo method which we 
use for their numerical solution. 



D. Classical model 

We want to stress that our perturbative theory of 
nonlinear coherent backscattering is not only valid for 
an atomic medium, but can be adapted to other kinds 
of nonlinear scatterers. In particular, the effect of in- 
terference between three amplitudes is always present 
in the perturbative regime of a small x^ 3 ^ nonlinearity. 
Specifically, we have also examined the following model: 
a collection of classical isotropic scatterers, situated at 
positions rj, i = 1, . . . , N. In analogy to the atomic 
model, we assume that the field scattered elastically by 
an individual scatterer at position ri is proportional to 
Ei/(1 + s\Ei\ 2 ), where Ei is the local field at ri, and s 
measures the strength of the nonlinearity. Writing Ei as 
a sum of the incident field, and the field radiated by all 
other scatterers, we obtain the following set of nonlinear 
equations: 



E, 



E; 



^ kr l3 i + s\E^' 



(60) 



Employing diagrammatic theory similar to the one out- 
lined above, we have checked that, in the ensemble aver- 
age over the positions rj, this model indeed reproduces 
the elastic components of the backscattering signal of the 
atomic model. We have checked that the results obtained 
from direct numerical solutions of the field equations l|60|) 
- averaged over a sufficiently large sample of single real- 
izations - agree with our theoretical predictions, in the 
perturbative regime of small nonlinearity s. In particular, 
the diagrams (d,e,f) of Fig. describing the interference 
contributions from the nonlinear propagation, are essen- 
tial to give the correct results. A more detailed analysis 
will be presented elsewhere. 



Furthermore, it remains to be clarified whether the 
diagrams (d,e,f) are also relevant for the description of 
propagation in homogeneous nonlinear media, into which 
linear scatterers are embedded at random positions. First 
studies of the resulting CBS cone have been presented 
in |20l I2II I , without taking into account interference be- 
tween three amplitudes, however. Experimentally, this 
question can be resolved by measuring the value of the 
backscattering enhancement factor 77: whereas 77 is basi- 
cally unaffected by the nonlinearity according to [2(1 l2ll ] 
(i.e. 7/ = 2 apart from single scattering), our equations 
(|44I46|I . with s proportional to the incoming intensity and 
to the coefficient of the nonlinear Kerr medium, pre- 
dict a significant change of 77 when varying the incoming 
intensity. 



IV. RESULTS 

We return to the atomic model, concentrating on the 
case of a slab geometry in the following. Using the equa- 
tions derived in Sec. IIII AIIIl"Cl we are able to calculate 
the backscattered intensity up to first order in the satu- 
ration parameter s. In this section, we will examine its 
dependence on the optical thickness b and detuning 5, 
for the scalar and vectorial case. The main quantity of 
interest is the backscattering enhancement factor rj. It is 
defined as the ratio between the total detection signal in 
exact backscattering direction divided by the background 
component. If we perform an expansion up to first order 
in s, we obtain 



V = 



L + C 
L 



71 (i) + ( ?? (i)_l)( 7c _ 7L ) s . 



(61) 



Here, rj 1 - ^ = 1 + /L^ is the enhancement factor in 
the linear case (i.e. the limit of vanishing saturation). 
If single scattering is excluded (e.g. in the h || h chan- 
nel), we have = 2. Increasing saturation changes the 
enhancement factor, and the present approach allows us 
to calculate the slope drj/ds of this change at s = 0. It 
is given by the difference between the nonlinear crossed 
and ladder contribution, normalized as follows: 



1L 



7c 



L - 

sLW ' 
C - C« 



(62) 
(63) 



Obviously, an important question is the domain of va- 
lidity of the linear expansion, Eq. (|61|l . Strictly speaking, 
this question can only be answered if we know higher or- 
ders of s. However, a rough quantitative estimation can 
be given as follows: if p\ (resp. P2+) denotes the prob- 
ability for a backscattered photon to undergo one (resp. 
more than one) nonlinear scattering event, the pertur- 
bative condition reads P2+ *C Pi- If we assume that 
all scattering events have the same probability (propor- 
tional to s) to be nonlinear (thereby neglecting the inho- 
mogeneity of the local intensity), we obtain p\ ~ (N)s, 
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FIG. 9: (color online). Normalized inelastic ladder and 
crossed contributions J^q, cf. Eqs. H62I63L for optical thick- 
ness b = 0.5, as a function of the laser detuning 6. Solid 
lines: polarized case (h\\h channel). Dotted lines: scalar case. 
For comparison, the corresponding elastic contributions (in- 
dependent of 8) are 7^ cl - ) = —7.04 (h \\ h), —6.53 (scalar), and 
7^ cl) = -9.56 (ft || ft), -18.8 (scalar). 



and P2+ — (N 2 )s 2 , where N denotes the total number of 
scattering events, and (...) the statistical average over all 
backscattering paths. Evidently, N and TV 2 are expected 
to increase when increasing the optical thickness b. For a 
slab geometry, we have found numerically that (N) oc b 
and (-/V 2 ) ex b 3 (in the limit of large b), concluding that 
the perturbative treatment is valid if sb 2 <C 1. Let us 
note that a similar condition also ensures the stability of 
speckle fluctuations in a nonlinear medium |43| . 

In Fig. El we show the inelastic ladder and crossed con- 
tributions 7^ m ' ) and for a slab of optical thickness 
b = 0.5 as a function of the detuning, 8 = lj — oj at , for 
the polarized (h \\ h) and scalar case. Since the opti- 
cal thickness is kept constant, the elastic quantities are 
independent of the detuning, and only the inelastic com- 
ponents are affected by 5, via the shape of the power 
spectrum P{lo') of the inelastically scattered light, see 
Eq. (|26p. The latter exhibits two peaks of width T, one 
of which is centered around the atomic resonance. The 
increase of the ladder term as a function of 8, which is 
observed in Fig.|!Ja) is due to initially detuned photons, 
i.e. uj = uj at + S, which are set to resonance (a/ ~ w a t) 
by the nonlinear scattering process. For these photons, 
the scattering cross section increases, which increases the 
contribution to the backscattering signal in the sum over 
all scattering orders - especially in the h || h case where 
single scattering is filtered out. The same effect also ap- 
plies for the crossed term, Fig. OJb), but here the de- 
phasing between the reversed paths due to the frequency 
change - which is more effective for higher values of the 
detuning - is dominant, leading in total to a decrease of 
7^ as a function of 8. The small ripples in Fig. |SJa), 
for the polarized case (solid line) at large 8, are due to 
numerical noise in the Monte-Carlo integration. 




FIG. 10: (color online). Normalized inelastic and elastic lad- 
der and crossed 7^ el (^ n ', cf. Eqs. 1621631 . for vanishing detun- 
ing 8 = 0, as a function of the optical thickness b = 0.5. Solid 
lines: polarized case (h\\h channel). Dotted lines: scalar case. 



Fig. shows the elastic and inelastic ladder and 
crossed contributions, as a function of the optical thick- 
ness, at detuning 8 = 0. The main purpose of this figure 
is to show the increase of the nonlinear contributions as 
a function of b, which is important to understand the do- 
main of validity of the present approach. The origin of 
this increase is simple to understand: for larger values of 
the optical thickness, the average number of scattering 
events increases, and so does also the probability that 
at least one of them is a nonlinear one. Thus, for an 
optically thick medium, even a very small initial satu- 
ration may lead to a large inelastic component of the 
backscattered light. Note, however, that the elastic and 
inelastic ladder contributions, Fig. llUf a.c). tend to cancel 
each other, such that their sum depends less strongly on 
b. Physically, this fact is related to energy conservation. 
The latter ensures that the total nonlinear scattered in- 
tensity - integrated over all final directions - vanishes even 
exactly, since the total outgoing intensity must equal the 
incident intensity (meaning a purely linear relationship 
between outgoing and incident intensity). 

Furthermore, we note that both the elastic and inelas- 
tic ladder components increase significantly slower in the 
polarized than in the scalar case (solid vs. dashed line). 
This is due to the fact that, as discussed in Scc. lIII 01 po- 
larization effects diminish the weight of nonlinear scatter- 
ing by approximately 2/3. Concerning the crossed com- 
ponents, Fig. llUf b.d). the difference is even stronger, due 
to the additional polarization-induced loss of contrast by 
a factor 3/4, in average. Please note that the vertical 
scale for the elastic crossed case, Fig. HUf d). is two times 
larger than in the other three cases: this reflects the effect 
of interference between three amplitudes, which renders 
the crossed component up to two times larger than the 
ladder. Concerning the inelastic component, Fig. llUf b). 
this effect is diminished by decoherence due to the fre- 
quency change at inelastic scattering. Here, crossed and 
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FIG. 11: (color online). Slope of backscattering enhancement 
factor, for the parameters of Fig. |U] (b — 0.5, left half) and 
Fig. 1101 (5 = 0, right half). Solid lines: polarized case (h\\h 
channel). Dotted lines: scalar case. 



ladder component are of similar magnitude. 

In Fig. ^2 we show the slope of the backscattering 
enhancement factor, which follows via Eq. (|61|) from the 
data shown in Figs. 151 andlTUl Fig. lTTT b) again points out 
the importance of even small saturation in the case of an 
optically thick medium. For example, in the scalar case 
at b = 2, increasing the saturation from s = to s = 0.01 
decreases the enhancement factor from 1.73 (< 2 due to 
single scattering) to 1.55. For very large 6, we find a lin- 
ear decrease of the slope. At the same time, however, the 
allowed domain of s <C 1 /b 2 shrinks to zero quadratically. 
This allows the enhancement factor to remain a contin- 
uous function of s, even in the limit b — > oo, where its 
slope at s — diverges. In order to make more precise 
statements about the behavior in the limit b — > oo, how- 
ever, it is necessary to generalize our theory to the case 
of more than one nonlinear scattering event. 

On the left hand side, Fig. Illf a) depicts the depen- 
dence of the enhancement factor on detuning, for b = 0.5. 
As already discussed above, the decrease of r\ with in- 
creasing S originates from the form of the inelastic power 
spectrum, which results in a stronger dephasing between 
reversed paths for larger detuning. Thus, the modifica- 
tion of the enhancement factor with the detuning, keep- 
ing fixed the linear optical thickness, is a signature of the 
nonlinear atomic response and has been experimentally 
observed in ref. |4£j. Let us stress, however, that in the 
cases shown in Figs. 151 and lllf a) and small detuning, the 
inelastic component gives a positive contribution to the 
backscattering enhancement factor. Hence, the observed 
negative slope of 77 originates from the elastic component, 
where the nonlinear crossed term is up to two times larger 
than the ladder, but with negative sign, see Fig. HUT cd). 

In order to observe an enhancement factor larger than 
two - and thereby demonstrate clearly the effect of in- 
terference between three amplitudes - it is therefore nec- 
essary to filter out the elastic component. In principle, 
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FIG. 12: (color online). Spectral dependence of the enhance- 
ment factor, for detuning 5 = and optical thickness b = 0.5 
(dashed line), 1 (solid), and 2 (dotted), in the h\\ h channel 
(a) and the scalar case (b). The vertical dashed line displays 
the position of the elastic S peak, which must be filtered out 
in order to observe an enhancement factor larger than two. 
The inset shows the power spectrum of the backscattered light 
(background component), which is almost identical with the 
single-atom spectrum. 
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FIG. 13: (color online). Spectral dependence of the enhance- 
ment factor, for detuning 5 — V and optical thickness b = 0.5 
(dashed line), 1 (solid), and 2 (dotted), in the h\\h channel 
(a) and the scalar case (b). The vertical dashed line displays 
the position of the elastic 5 peak, which must be filtered out in 
order to observe an enhancement factor larger than two. The 
inset shows the power spectrum of the backscattered light 
(background component), revealing the amplification of the 
on-resonance peak with respect to the symmetric single-atom 
spectrum. 
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this can be achieved by means of a spectral filter, i.e. by 
detecting only photons with a certain frequency a/, dif- 
ferent from the laser frequency uj. Thereby, it is possible 
to measure the spectral dependence of the backscatter- 
ing enhancement factor, see Fig. ^21 Here, the upper (a) 
and lower (b) half depict the polarized (h || h) and scalar 
case, respectively, for vanishing laser detuning, 5 = 0. 
Evidently, the largest values of the enhancement factor 
are obtained if the final frequency approaches the initial 
one, since then the dephasing due to different frequencies 
vanishes. In the scalar case, the value of the enhance- 
ment factor in the limit 8' — > is completely determined 
by the relative weights between the one-, two- and three- 
amplitudes cases shown in Fig.EJ cf. Eq. As evident 
from the dashed line in Fig. I12f b) , already at the rather 
moderate value b = 0.5 of the optical thickness, the three- 
amplitudes case is sufficiently strong in order to increase 
the maximum enhancement factor above the linear bar- 
rier i] = 2. With increasing optical thickness (and, if 
necessary, decreasing saturation parameter, in order to 
stay in the domain of validity of the perturbative ap- 
proach, see above), the number of linear scattering events 
increases, which implies that the three-amplitudes case 
increasingly dominates, see Fig.^J In this limit, the en- 
hancement factor approaches the maximum value three. 
At the same time, however, a larger number of scattering 
events also leads to stronger dephasing due to different 
frequencies, uj' =/= uj. This results in a narrower shape of 
r\ as a function of uj' for larger optical thickness. Nev- 
ertheless, as evident from Fig. I12f b). the enhancement 
factor remains larger than two in a significant range of 
frequencies uj' . The same is true for the polarized case, 
Fig. I12f a). However, here the enhancement factor can- 
not exceed the value 2.5, due to the polarization- induced 
loss of contrast. At the same time, the optical thickness 
has less influence on the maximum enhancement factor 
at 8' = 0, since single scattering, Fig. 0Ja) - and partly 
also the two-amplitudes case, Fig.^Jb) - are filtered out, 
so that interference of three amplitudes already prevails 
at rather small values of the optical thickness. 

In Fig. E| the influence of an initial detuning (here: 
8 = r) is displayed. Basically, the above conclusions re- 
main almost equally valid for the detuned case. A small 
difference is seen in the scalar case, Fig. Ilijf b). where 
the maximum of ry(<5') is found slightly below 8. This is 
due to the fact that the weight of single scattering in- 
creases with increasing §'. Furthermore, the inset reveals 
that the power spectrum of the backscattered light dif- 
fers from the single-atom spectrum, Eq. I|26|l . where the 
two peaks at 8' = and 8' = 28 are equally strong. In 
the multiple scattering case, the on-resonance peak at 
8' = is amplified, since the scattering cross-section is 
larger for photons on resonance. As already mentioned 
above, see the discussion of Fig. [5] this increases the total 
contribution to the detection signal (in the sum over all 
scattering paths) - especially in the polarized case, where 
single scattering is filtered out. 



V. CONCLUSION 



In summary, we have presented a detailed diagram- 
matic calculation of coherent backscattering of light from 
a dilute medium composed of weakly saturated two-level 
atoms. Our theory applies in the perturbative two- 
photon scattering regime (s -C 1 and sb 2 <C 1), where at 
most one nonlinear scattering event occurs. The value of 
the backscattering enhancement factor is determined by 
the following three effects: firstly, due to the nonlinearity 
of the atom-photon interaction, there may be either two 
or three different amplitudes which interfere in backscat- 
tering direction. This implies a maximal enhancement 
factor between two and three for the nonlinear compo- 
nent, where the value three is approached for large optical 
thickness. However, since the contribution from nonlin- 
ear scattering has a negative sign, the total enhancement 
factor (linear plus nonlinear elastic and inelastic compo- 
nents) is reduced by the effect of three-amplitudes inter- 
ference. Only if the elastic component is filtered out, a 
value larger than two can be observed. 

Secondly, a loss of coherence is implied by the change 
of frequency due to inelastic scattering - like in the case 
of two atoms |32j . The random frequency change leads to 
different scattering phases - and hence on average deco- 
herence - between reversed paths. Finally, a further loss 
of contrast is induced by nonlinear polarization effects - 
even in the h \\ h channel, which exhibits ideal contrast 
in the linear case. Nevertheless, the enhancement factor 
remains larger than two in certain frequency windows of 
the inelastic backscattering signal. Thus, it is experimen- 
tally possible to clearly identify the effect of interference 
between three amplitudes - provided a sufficiently narrow 
spectral filter is at hand. 

A natural way to extend this work is to give up the 
perturbative assumption, and admit more than one non- 
linear scattering event. This is necessary in order to de- 
scribe media with large optical thickness, even at small 
saturation. Since the number of interfering amplitudes 
increases if more than two photons are connected by 
nonlinear scattering events, we expect the occurrence of 
even larger enhancement factors in the nonperturbative 
regime - especially in the case of scatterers with posi- 
tive nonlinearity, i.e. for scatterers whose cross section 
increases with increasing intensity. 

Furthermore, the relation between coherent backscat- 
tering and weak localization in the presence of nonlinear 
scattering remains to be explored. Does a large enhance- 
ment of coherent backscattering also imply a strong re- 
duction of nonlinear diffusive transport? If the answer is 
yes - as it is the case in the linear regime - this implies 
that wave localization can be facilitated by introducing 
appropriate nonlinearities. 
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APPENDIX A: MONTE-CARLO SIMULATION 



-TT( c )Cp f ? p ? w 

. ±± y^u: ^Vl ^Wl ^w? } ^ 



As discussed in Sec. MI CI the incorporation of polar- 
ization effects requires to take into account the projection 
of polarization vectors in the corresponding scalar equa- 
tions. For the inelastic ladder component, insertion of 
the polarization term, Eq. Q55p. into the scalar expres- 
sion, Eq. J2HJ), yields: 



0, if m = I = 0, 
2, ifn = m = 0, ^>0 
2, ifn = Z = 0, m>0, 
4, otherwise 

(A3) 



with the polarization vectors of the 'reversed' photons: 



= A 



v m ,v m+ i 



■ A VliV2 ££>, 



A Wl Wl+1 . . . A WljW2 e£. 



(A4) 



1=0 



The elastic nonlinear scattering components follow sim- 
ply by inserting — 2<5(u>'— ui) instead of the inelastic power 
spectrum P{u>') in the above Eqs. (|A1IA3|I . 

The nonlinear propagation term is obtained by insert- 
ing Eq. iJSSJl into Eq. gljl: 



du^ . . . du n e 



rfvi . . . dv m e 



-ui, z /e 



j [ <Sa;G w (Ui, Ui+i) 



\i=l 



L 



(2, prop) 
il 



j=l J \i=l ) m=l ^ V 



/ dw! ...dwi e -"*-/* f TJ |5^G^(w k ,w k+1 ) 



fe=l 

1 if n = m = 0, 



E 



/'k^ejxr (Ai) 

2 2 it n > or m > 0, 



with u n+ i = v m+ i = wi + i = r. Furthermore, the po- 
larization vectors are given by: 



€ v — ^v m .v m _|_i • • ■ ^vi,V2^L? 
E w — ^. WljW j, 1 . . . A Wl;W2 6x). 



n^He.e.e^x J if m _ Z - 0, 
2 otherwise, 



Here, the nonlinear event takes place between ui 
and ui + i. Correspondingly, J Ul+1 denotes the one- 
dimensional integral on a straight line between these 
points, and Uo = Ui — Ui, z e z and u n+ i = u n — u n ^ z e z are 
defined as the points where the photon enters or leaves 
the medium, respectively. The three cases Fig. EIa,b,c) 
correspond to < I < n, I — n, and I — 0, respectively. 
The polarization vectors t\ and €3 participating in the 
nonlinear event, cf. Fig. |SJ are obtained as: 



(A2) 



^ui+i,ui ■ • • A U2 Ul e/ J , 

^Ul,Ul+l 



e 3 — A UliUl+1 . . . A Un l jUn ei3. 



(A6) 
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Finally, to obtain the interference component C^' plop \ 
the last term in Eq. I|A5|1 must be replaced by: 

'2 if n> l,(m,0 + (0,0), 
n (i ' prop )(ei,e^,e 3 ) x \ 1 ifn>l,m=l = 0, + 
otherwise 

4 if Z = 0,n > l,m > 0, 
2 if / = 0,n = l,m > 0, 
2 if < Z < n, 
otherwise, 

(A7) 



with 



£3 



i ui,ui+l • ■ 



(A8) 



The first term, n^ i,prop \ equals the ladder component 
minus single scattering (n = 1), whereas the second 
one, n( c,prop ), describes the additional crossed diagrams 
shown in Fig. EJd-f). Here, the case < I < n corre- 
sponds to Fig. Efd), where the nonlinearity occurs be- 
tween two scattering events. The remaining diagrams, 
Fig-IUe,f), correspond to I = 0. Here, the case m = 
('pump photon from the coherent mode') does not con- 
tribute, since then the diagrams Fig. HJe,f) are identical 



to Fig. HJb,c). Furthermore, if n = 1 ('probe photon 
singly scattered'), the two diagrams Fig.HJe) and (f) be- 
come identical. In this case, we obtain a factor 2, whereas 
the sum of diagram (e) plus diagram (f ) yields 2 + 2 = 4 
in the case n > 1. 



Numerically, we solve the above integrals by a Monte- 
Carlo method. Here, we proceed as follows: for 
Eqs. (|A1IA3(1 . at first position r and frequency ujd of the 
inelastic scattering event are chosen randomly. Start- 
ing from r, three photons are launched, two with fre- 
quency lul and one with frequency ljd- After each scat- 
tering event, the length r of the next propagation step 
is determined randomly according to the distribution 
P(r) — exp(— r /£)/£, whereas the direction is chosen uni- 
formly. After all photons have left the medium, the triple 
sum over n, m, and / is performed, taking into account 
the projection of the polarization vectors. For the nonlin- 
ear propagation term, Eq. I|A5|) , at first the probe photon 
(path: Ui,...,u n ) is propagated, starting in the laser 
mode 1cl,£l. Then, the pump photon is launched from 
a randomly chosen position v m on the path of the probe 
photon. Finally, the projection of polarization vectors is 
performed separately for each given path. 
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